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Abstract. We present deep V,I photometry of two 9.4' x 9.4' field in the outer regions of the M 33 galaxy. We obtain a robust 
detection of the luminosity of the Red Giant Branch Tip (J TRGB = 20.72 ± 0.08) from which we derived a new estimate of 
the distance modulus of M 33, (m - M) = 24.64 ± 0.15, corresponding to a distance D = 847 ± 60 Kpc. By comparison of 
the color and magnitude of the observed Red Giant Branch stars with ridge lines of template globular clusters we obtained 
the photometric metallicity distribution of the considered fields in three different metallicity scales. The derived metallicity 
distributions are very similar over a range of distances from the galactic center 10' < R < 33', and are characterized by a well 
defined peak at [M/H] - -0.7 ([Fe/H] - 1.0, in the Zinn & West scale) and a weak metal-poor tail reaching [M/H] i - 2.0. 
Our observations demonstrate that Red Giant Branch and Asymptotic Giant Branch stars have a radial distribution that is much 
more extended than the young MS stars associated with the star-forming disc. 
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1. Introduction 

1 The accurate determination of the distance to Local Group 
galaxies is critical for the establishment of a reliable extra- 
galactic distance scale. In particular, local late-type galaxies as 
the Large Magellanic Cloud, M 31 and M 33 are the sites of 
choice for the calibration of several secondary distance indica- 
tors. 

The Sc II-III spiral M 33 (NGC 598 or Triangulum 
, galaxy) is the t hird-brightest member of the Local Group (see 
IVan den Berghl Ho00, for a review). In spite of that, while 
Cepheid variables wer e discovered as early as 1920 in this 
galaxy (Hubble, 1926), reliable d istance estimates have be- 
come available only after 19 80 Jsandage & Carlson! 1 19831 
IChristian & Schommerl Il987l) and differences of the order of 
~ 0.3 mag can be found even amon g the most recent estimates 
of the M 33 distance modulus (see McConnachie et al 1 12004 
and references therein). 

Concerning the stellar populations of M 33, most of the 
studies have centered their attention on the massive (young) 
stars popula t ing the disc (s e e, for example IVan den B ereh, 
Il99ll l2000t lUrbaneia et all l2002l and ref erences therein) 
and/or on the innermost regions of the galaxy jMighell & Richl 
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ll995llStephens & Froge3.l2002l) . On the other hand, only a few 
authors have provided some insight into the outer regions and 
the (possible) halo population ( see JMould & KristianI (l986), 
IWilson et all J 199(1 iDayidg el d2003l) and the HST-WFPC2 
study bv lKim et al.N2002l) ). 

Here we present V and I photometry (reaching V ~ 24.5) 
of two fields located at ~ 15' and ~ 28' from the center of 
M 33, to the North- West of the main body of the galaxy, ap- 
proximately in the direction of the minor axis (see Fig. 1). The 
Tip of the Red Giant Branch (TRGB) is cleanly detected with 
our data and we use it to derive a robust estimate of the dis- 
tance modulus of M 33 adoptin g the RR Lyrae-ind ependent 
calibration prov i ded by iBellazzini. Ferraro & Pancind J200ll) . 
Bellazzin i et al.l d2004l) . This is the main aim of the present 
study which is part of a large programme devoted to the de- 
termination of homogen eous distances for most of the galaxies 
of the Local Group (see Bellazzi ni et all 12002). The metallic- 
ity distribution of the RGB stars in the considered fields is also 
obtained from photom etric estimates, following the method by 
IBellazzini et al.U2003l) . Shortly after this p aper was submitted, 
a preprint was posted JTiede et all l2004 hereafter T04) pre- 
senting the analysis of a 6.8' x 6.8' field located in the South- 
Eastern region of M 33 at ~ 20' from its center. The analysis is 
very similar to that performed here and the results are in excel- 
lent agreement, as we will show below. 

The structure of the paper is as follows: in Sect. 2 we de- 
scribe the observational material, the data reduction process 
and the photometric calibration; in Sect. 3 we present the color- 
magnitude diagrams (CDMs), our TRGB estimate of the dis- 
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Fig. 1. The position of the observed fields over-plotted on a 
DSS-II image of M33. The dimension of the image is 40' x 33'. 



tance modulus, and the (photometric) metallicity distributions. 
The stellar content of the considered fields is also briefly dis- 
cussed. A brief summary and the main conclusions are reported 
in Sect. 4. 

2. Observations and data reduction 

2.1. Observations 

The observations were carried out on September 7 and 8, 2002, 
with the imager/spectrograph DoLoRes at the 3.52m Italian 
telescope Telescopio Nazionale Galileo (TNG) in La Palma 
(Canary Island, Spain). DoLoRes is equipped with a 2048 x 
2048 pixel thinned and back-illuminated Loral CCD array with 
a total field of view of 9'. 4 x 9'. 4, and a scale of 0.275 "/px. 
The seeing was remarkably good and stable, ranging between 
0.8" and 1.1" FWHM for all the scientific images. 

Two fields were imaged along the minor axes in V and I 
filters. The locations of the fields are shown in Figure 1 and 
listed in Table 1 along with the date, the filter and the expo- 
sure time of all the scientific exposures. A total of 7 V and 6 I 
600 s exposures have been acquired for the inner field (Field 1, 
hereafter Fl) and 8 V and 6 I 600 s for the outer field (Field 2, 
hereafter F2). The distance of the field center from the center 
of the galaxy is ^ 15' for Fl and 28' for F2. 

2.2. Data analysis 

All the images were corrected for bias and flat-field using stan- 
dard IRAF procedures. The relative photometry was carrie d 
out with the PSF-fitting code DoPhot JSchechter et alJ dl993h . 
Each frame was independently reduced. A 3 <x threshold above 
the background noise was adopted for the search for sources 
on the frames and the spatial variations of the PSF were mod- 
eled with a quadratic polynomial function. Only the sources 
classified as stars by the code were retained. All the V and 
I catalogues of a given field were reported to the (instrumen- 
tal) photometric system of the best-seeing image acquired un- 
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Table 1. Log of the observations 



Field RA (J2000) 


Dec (J2000) 


Date 


Filt. 


Exp T(s) 


Fl 01 32 44.85 


+30 45 02.5 


Sept,7,02 


V 


3x600 






Sept,7,02 


1 


3x600 






Sept,8,02 


V 


4x600 






Sept,8,02 


I 


3x600 


F2 01 31 44.97 


+30 45 00.0 


Sept,7,02 


V 


5x600 






Sept,7,02 


1 


4x600 






Sept,8,02 


V 


3x600 






Sept,8,02 


I 


2x600 



der photometric conditions. The catalogues were then cross- 
correlated, the magnitudes were averaged and the standard de- 
viation adopted as the photometric uncertainty of the individ- 
ual stars. In the final catalogue of instrumental magnitudes and 
positions we have retained only the sources that have at least 
three valid measures of the magnitude for each passband (e.g. 
at least 3 V and 3 I). Moreover, all the stars with an associated 
error (either in V or I magnitude) larger than three times the av- 
erage uncertainty at their magnitude level were excluded from 
the sample. The final Fl and F2 catalogues contains 26399 and 
1937 sources, respectively. Accurate aperture corrections have 
been obtained for each field on a few tens of bright and isolated 
stars. 

The stellar crowding is quite low in the considered fields 
that sample external low-surface brightness regions of the 
galaxy. Even in Fl the average stellar density is as low as 
0.08 stars /circsec 2 (considering all the detected stars). For the 
applications presented in this paper the effects of incomplete- 
ness are not an issue since we always compare subsets of stars 
that are homogeneous in magnitude and color (see Sect. 3.3 
and 3.4 below). No attempt is made to determine the true 
Luminosity Functions of the observed sequences. The TRGB 
level, the main target of the present study, occurs more than 3 
magnitudes above the limiting magnitude, a range in which the 
considered sample is likely ~ 100% complete. Since the pho- 
tometric uncertainties of individual stars (in each passband) are 
empirically estimated as the standard deviations of 3 to 8 inde- 
pendent measures of the magnitudes there is no need for artifi- 
cial star experiments to characterize photometric errors. 

2.3. Photometric calibration 

The absolute calibratio n has been ob tained from several re- 
peated observations of lLand olt (1992) standard fields, includ- 
ing all the st ars listed in th e extended catalog of calibrators 
provided by Stetson (2000). The coefficients of atmospheric 
extinction (C ext ) were directly obtained by repeated observa- 
tions of the same standard field at different air mass. The final 
calibrating equations and C ext are shown in the upper panels of 
Figure 2. 

To check the accuracy (and reproducibility) of the photo- 
metric calibration we reduced one V and one I 30 s frames 
centered on the globular cluster NGC 6779, acquired during 
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Photometric Calibration 




Fig. 2. Upper panels: difference between tabulated magnitudes 
(V,I) and instrumental magnitudes ( v,i) vs instrumental color 
index (v-i) for the observedl Landoltl (H"992) standard stars. The 
calibrating relations are plotted (solid lines) and the corre- 
sponding equations are reported along with the RMS of the 
linear fits and the value of extintions coefficients. Lower pan- 
els: comparison of the photometry of NGC 6779, obtained from 
data acquired during the same run of the M 33 data and cali- 
brated with th e relations abov e, and the independent photome- 
try bvlRosenber g et alJ (E"000). 



the second night of the run. The catalogue was calibrated 
with the relations shown in Fig. 2 and the final photometry 
was compared with the photometry of the cluster provided by 
iRosenberg et al.l lEoOO). The results of the comparison (shown 
in the lower panels of Fig. 2) demonst rate that our pho t ometr y 
is in excellent agreement with that of Rosenberg et all d2000l) . 
an the characteristic uncertainty of our absolute calibration is 
of the order of < +0.02 mag. 

3. Results 

3.1. Color-Magnitude Diagram 

The Color Magnitude Diagrams (CMD) of the observed field 
are presented in Fig. 3. The CMD of Fl is dominated by a 
wide RGB sequence, running from / ~ 20.6 down to the 
limiting magnitude of the photometry. The exact location of 
the TRGB (as derived in Sect. 3.2, below) is reported in the 
CMDs as a horizontal line to the red of the observed RGB. 
The sources brighter and redder than the RGB Tip are likely 
bright Asymptotic Giant Branch (AGB) stars, probably asso- 
ciate d with an interm ediate-age (and/or metal-rich) population 
(see lDavidgell2003l for an extensive discussion). The presence 
of a conspicuous plume of intermediate-young Main Sequence 



Fig. 3. Color Magnitude Diagrams of the two observed fields. 
The errobars show the average photometric uncertainty as a 
function of magnitude; the horizontal lines marks the position 
of the TRGB as determined in Sect. 3.2. 



(MS) stars reaching / ^ 19.4 is also evident on the blue side of 
the CMD, with a sharp blue edge at V - I =* 0.0. In the CMD 
of F2 only the RGB and AGB population are discernible, while 
the blue MS stars are completely absent. Star counts on the 
RGB indicate that the stellar density drops by a factor > 20 
from Fl to F2. 

The contamination by foreground Galactic stars is negligi- 
ble for the purpose s of the present study. The Galactic model by 
lRobinetalJ (l2003) predicts fewer than 200 Galactic stars in the 
observed fields within the color and magnitude ranges spanned 
by the CMDs of Fig. 3. By inspectio n of the synthetic CMD 
obtained with the Robi n et al.l J2003) model we conclude that 
even in the poorly populated F2 the only region of the CMD 
that may be significantly affected by Galactic contamination is 
that which hosts bright AGB stars, e.g. V-I > 1.8 and V i 20.6 
(see Sect. 3.4, below). 

3.2. TRGB distance 

The use of Tip of the Red Giant Branch (TRGB) as a 
standard candle is now a widely used technique to esti- 
mate flie_d^stanc£_Jo__gakdxies of an y morpholo gical type 
(see |Lee . Freedman & Madora Il993t iMadore & Freed man. 
Il995l Il998l IWalkerl 120031 for a detailed description of 
the method, recent reviews and app lications). The underly- 
ing physics is well und e rstood JMadore & Freedmarl [l998; 
ISalaris. Cassisi & Weissl l2002h and t he observational pro- 
cedure is operationally well defined JMadore & Freedmanl 
1995). The key observable is the sharp cut-off occur- 
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ring at the bright end of the RGB Luminosity Function 
(LF) that can be easily detected wit h the application of 
an edge-detector filter (Sobel filter. [ Madore & Freedmanl 
1 19951 ISakai. Madore & FreedmanL 1 19961) or b y other (gener - 
ally parametric) technique s (see, for example iMendez et all 
120021 iMcConnachie et all I2004I) . The necessary condition 
for a safe application of the technique is that the observed 
RGB Luminosity Function should be well populated, with 
more than ~ 100 stars within 1 mag from t he TRGB 
(Madore & Freedman[lr99llBellazzini et allEool). 

The F2 sample is not sufficiently populated for a safe ap- 
plication of the method while the Fl sample clearly fulfils the 
above criterion (there are more than 2500 RGB stars within 1 
mag from the TRGB), hence we limit the TRGB research to 
Fl. As a first step, to limit the range of metallicity of the stars 
involved in the TRGB detection, we select RGB stars by colo r 
following the approach adopted bv lMcConnachie et alJ (2004). 
The adopted selection includes the main bulk of the RGB pop- 
ulation and it is shown in the lower left panel of Fig. 4. The 
logarithmic LF is presented as an ordinary histogram and as 
a generalized histogram (e.g. the histogram convolved with a 
Gaussian with standard deviation equal to the photometric error 
at the given magnitude, see lLaird et allll988HBellazzini et all 
l2002l for definitions and references) in the upper left and up- 
per right panels of Fig. 4, respectively. The sharp cut-off is 
an obvious feature of both representations of the LF and is 
easily detected by the Sobel filter (Fig. 4, lower right panel). 
As usual, the peak of the filter response is taken as the best 
estimate of the TRGB location and the Half Width at Half 
Maximum of the same peak is taken as the associated uncer- 
tainty, I TRCB = 20.72 ± 0.08. If we consider the most recent 
estimates in the literature, our valu e is ~ 2.2cr larger than that 
found bv M cConnachie et alJd2004l I TRGB = 20.54 ± 0.01; but 
these authors provide only a formal error on their est i mate), 
and ~ 1 — 2cr lower than the estimates bv iKim et alJ (|2002, 
jTrgb _ 20.82 - 20.92 ± 0.05 depending on the considered 
field), e.g. it is bracketed by the two quoted results. On the other 
hand our estimate is in excellent agreement with that obtained 
by T04 (I TRCB = 20.75 ± 0.02). 

We a dopt E(B — V) = 0.04, ac cording to the reddening 
maps by Schlegel etafl dl998l) and iBurstein & Heilesl Jl 984) 
and A j = 1 J6E(B - V), according to lDean. Warren & Cousins! 
( 1978). We note however that most of the other available esti- 
mates of the foreground redden ing cluster around E(B - V) = 
0.08 (see lVan den Berghlll99ll and references therein). To ac- 
count for this, we report also the results we obtain adopting 
E(B -V) = 0.08 (see Table 2, below). Note that the effects of 
this different assumption are small either on the final distance 
modulus (e.g. 0.05 mag) or the average metallicity (< 0.15 

dex; see Tab. 2) . According to the detailed dust maps of M 33 

I ' ■ I 1 1 

by Hippelein et al. (2003) the effect of the intrinsic extinction 

should be negligible in the fields con sidered here. 

In lBellazzini. Ferraro & Pancinol ll200ll) we have provided 
a robust zero-point to the calibrating relation providing the ab- 
solute I magnitude of the tip (Mj RGB ) a s a function of metallic- 
ity ([Fe/H], in the IZinn & Wesil d 1984 scale, hereafter ZW), 
based on the geometric dis t ance t o the cluster u> Centauri ob- 
tained by Thomp son et all {2001 ) using the double-lined de- 




Fig. 4. Detection of the TRGB. The CMD in the lower left 
panel shows the adopted selection, e.g. the stars enclosed by the 
two diagonal lines. The arrow marks the position of the TRGB. 
The upper panels display the logarithmic LF of the upper RGB 
as an ordinary histogram (left) and as a generalized histogram 
(right). The thick lines marks the position of the TRGB, the 
thin lines enclose the associate uncertainty range. Lower panel: 
Response of the Sobel's filter to the observed LF. 

tached eclipsing binary OGLE-17. This calibration is fully in- 
dependent of the usual (Cepheid b ased an d/or RR Lyrae based) 
distance scales. In Bell azzini et al.lll2004l) we have extended the 
calibration to Near Infrared passbands and we refined our I cal- 
ibration providing also the relation for Mj RGB as a function of 
the glo bal metallicity (I'M/HI, see ISalaris. Chieffi & Straniero 
(1993) and lFerraro et al 1 dl999i) for definitions and discussion) 
that we adopt in the present analysis. 

Since the distance modulus derived from l TRGB is weakly 
dependent on metallicity, and our metallicity estimates (ob- 
tained by comparison with template RGB ridge lines, see be- 
low) depend on the assumed distance modulus, we adopted 
an iterative method to find simultaneously the two quantities 
searched for. First we derived a preliminary distance modulus 
adopting Mj RCB = -4.04, then we derived a median metallic- 
ity of the considered population as described in Sect. 3.3 below, 
and we obtained a refined estimate of the modulus using the 
obtained median metallicity (\M/ H] _ -0.75, see below) a s an 
input for the calibrating relation bv lBellazzini et al] d2004l) : 



M 



TRGB 



= 0.258[M///] 2 + 0.676[M///]- 3.629 (+0.12). (1) 



The process converged to stable values of the distance mod- 
ulus and of the median metallicity in 2-3 iterations, indepen- 
dently of the assumed reddening and/or the considered metal- 
licity scale. Our final estimate (for E(B - V) = 0.04 and 
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[M/H] med = -0.75) is (m - M)o = 24.64 ± 0.15, where all 
the sources of uncertainty have been taken into account. The 
corresponding distance estimate is D — 847 ± 60 Kpc. 

3.3. Metallicity 

We derive the Metallicity Distributions (MD) of the studied 
fields from the Color and Magnitude distribution of RGB stars 
(transformed to the absolute (V - 7)o vs. M/ plane adopting 
the reddening and distance modulus described above) by in- 
terpolation on a grid of RGB ridgeli nes of template globula r 
clusters, adopting the same scheme as Bellazzi ni et al.l (12003). 
Essentially the same approach is adopted also by T04. The 
MDs are obtaine d in different metallicity s cales, e.g. the ZW 
scale, the scale bv lCarretta & Grat ton ( 1997, hereafter C G) an d 
the global metallicity scale described in Fer raro et alJ ( 1999), 
to make the comparison with other studies easier. While the 
individual photometric metallicities provided by the adopted 
procedure may be quite uncertain, the overall metallicity dis- 
tribution and its average properties are sufficiently well char- 
acterized to provide interesting insights into the stellar popu- 
lation under analysis and it has been widel y used, in recent 
years , in the study of resolved galaxies (see iBellazzini et all 
2003, for details, discussion and references). The color distri- 
bution of RGB stars should depend - to a lesser extent - also on 
the age distribution of the underlying population. T04 studied 
this problem in the case of M 33 by mean of synthetic CMDs 
drawn from theoretical evolutionary tracks and concluded that 
"...the ages of the RGB stars are not likely to significantly affect 
the derived MDs...". They estimate that the maximum expected 
shift of the peak of the MDs is of the order of 0.1-0.2 dex. 
Finally, the young MS population observed in Fl - and not in 
F2 - cannot affect the comparison between the MDs of the two 
fields since their evolved counterparts should be negligible in 
number and do not fall in the selection box we adopt to obtain 
MDs (see Fig. 5 and 7, below). When dealing with photometric 
metallicities it should be kept in mind that the underlying age 
distribution may affect the derived MDs. 

In Fig. 5 a direct comparison between the observed RGBs 
of Fl and F2 and the adopted template ridge lines is presented. 
It is immediately clear that the large majority of M 33 RGB 
stars (in both fields) are enclosed within the ridge lines of 
M 5 ([Fe/H] cc = -1.11; [M/H] = -0.90) and of 47 Tuc 
([Fe/H] cc = -0.70; [M/H] = -0.60). In the CMD of F2 an 
anomalous clustering of stars can be noted around the ridge 
line of NGC 6553 ([Fe/H] cc = -0.16; [M/H] = -0.06) that 
has no counterpart in the CMD of Fl. However, as can be ap- 
preciated from Fig. 3, such red stars are just above the limiting 
magnitude at their color. For this reason we will not discuss 
in detail this feature in the following. Deeper photometry is 
needed to firmly assess the possible presence of an excess of 
very red RGB stars in this region. 

In Fig. 6 we present the MDs (in the form of generalized 
histograms) as a function of the global metallicity for Fl (upper 
left panel) and F2 (upper right panel). To study in finer detail 
the radial behavior of he MD we split Fl in two subregions (A 
and B) of similar area. In particular Fl A contains all the Fl 



Fl F2 




(v-i)„ (v-i) 



Fig. 5. The CMDs of the RGB stars in Fl (left panel) and F2 
(right panel) are compared to the grid of RGB ridge lines of 
template globular clusters we adopted to derive the metallicity 
distributions. From blue to red, the template clusters are: M 92, 
M 13, M 5, 47 Tuc, NGC 6553 and NGC 6528. The thick box 
encloses the RGB stars actually selec ted for the deri vation of 
the metallicity distribution (see IBellazzini et all 120031 for de- 
tails). 



stars less distant than 15' from the center of M 33, while the 
Fl stars with r > 15' are assigned to Fl B. The MDs of Fl A 
and Fl B are plotted in the lower left and lower right panels 
of Fig. 6, respectively. Note that if we exclude from the inter- 
polation the stars with M/ > -2.5, i.e. in the range where the 
sensitivity of color to metallicity is lower, the obtained MDs 
are unchanged. This experiment also demonstrates that the de- 
rived MD are not sensitive to the effect of incompleteness, as 
expected. The average properties of all the considered MDs are 
summarized in Table 2, for two different assumptions of E(B- 
V). There are a number of considerations emerging from the 
inspection of Fig. 6 and Table 2: 

1. The MD of all the considered fields shows a strong peak 
at [M/H] =i -0.7. This justifies our assumption of the 
median metallicity as the characteristic value of the domi- 
nant population in our determination of the TRGB distance. 
A sparsely populated tail of metal-poor stars extending to 
[M/H] < -2.0 is also present in all the presented MDs. 
This general similarity over large areas of the galaxy (Fl 
covers a range of galactocentric distances from ~ 2.4 Kpc 
to ~ 5 Kpc, F2 from ~ 5.5 Kpc to ~ 8.2 K pc) is in agree- 
ment with the results bv lKimetalJ( 2002) and is reminis- 
cent of what is observed in M 31 (IBellazzini et all Eool . 
Our MDs are very similar to those obtained by T04. 
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-2 -1 -2 -1 

[M/H] [M/H] 



Fig. 6. Metallicity Distributions (continuous lines) of the Fl 
field as a whole (upper left panel), of F2 (upper right panel), 
of the inner part of Fl (10' < R < 15'; Fl A, lower left panel) 
and of the outer part of the same field (15' < R < 20'; Fl B, 
lower right panel). The dotted line in the upper left panel shows 
the instrumental response of the method, the MD that would be 
obtained if the width of the RGB was entirely due to the pho- 
tometric errors. 

2. The bell-shaped curve plotted as a dotted line in the up- 
per left panel of Fig. 6 displays the response of the adopted 
interpolation scheme to a Simple Stellar Population (SSP, 
i.e. a populatio n of stars having the same a ge and chemical 
composition, iRenzini & Fusi Peccii Il988i) observed under 
the same conditions as our real data. It has been obtained 
adopting, as input for the interpolation, a synthetic RGB 
population whose color width is entirely due to the pho- 
tometric errors. The I magnitude of the "synthetic" stars 
is extracted from the observed RGB luminosity function, 
the V-I color is obtained from the average ridge line of the 
observed RGB plus a photometric error drawn at random 
from a Gaussian distribution having cr equal to the average 
observed photo metric uncertainty at the considered magni- 
tude (as done in Bella zzini et a ll 120021) . 

To obtain the true width of the underlying MD the de- 
scribed "instrumental response" should be deconvolved 
from the observed MD (continuous line). The main peak 
of the observed MD is well fitted by a Gaussian distribu- 
tion with cr = 0.34 dex, while the instrumental response 
curve is well approximated by a Gaussian distribution with 
cr = 0.23 dex. It may be concluded that the true intrinsic 
dispersion of the main peak of the MD is cr ^ 0.25 dex. 

3. By comparison with ridge lin es of template globular clus- 
ters, iMsinil^Kriliianl estimated < [Fe/H] Z w >- 
-2.2 ± 0.8 for a field at a galactocentric distance sim- 



ilar to F2. On the other hand, we find, for both Fl 
and F2, < [Fe/H] zw >^ -1.03 + 0.40 in ex - 
cellent agreement with the re s ults b y iDavidg d (|2003), 
ICuillandre. Lequeux & Loinard £l999) and T04. The dif- 
ference is partly justified by the different assumptions about 
distance ((m-M ) = 24.8 instead of our (m-M ) = 24.64). 
However, even adopting their distance modulus we find 
< [Fe/H]zw >— -1. 22 + 0.40, much more me tal-rich than 
what was found by iMould & KrisTkrl (l986). We tenta- 
tively ascribe this difference to a possible problem in the ab- 
solute calibration of Mould & Kristian's photometry. This 
hypothesis is confirmed by the results of T04. 

4. Assum i ng the same dista nce modulus as lMould & Kristiai] 
ill 9861) . iKim eTail d2002l) obtain -0.61 << [Fe/H] zw >< 
-0.86 for 10 fields covering (approximately) the same ra- 
dial range than our Fl. The agreement with our results is 
much better than with Mould & Kristian's one, still the dif- 
ference is not negligible and it is not justified by the differ- 
ent distance modulus assumed (a larger distance modulus 
should im ply a brighter and hence more metal-poor RGB). 
However IKim et all d2002h derived their mean metallicity 
from the mean (V - 7)o color at M/ = -3.5. We consider 
our median/mean metallicities, derived from all the RGB 
stars brighter than Mj = -2.0 and based on an accurately 
checked photometric calibration, as more robust and safer 
than those bvlKim et all ilOO'A . 

5. lKimetalJ (2002) found a weak radial gradient in the mean 
metallicity (0.2 dex) in the range 1 < R < 5 Kpc (see also 
T04). We find no sign of variation of the mean metallicity 
in the range 2.4 < R < 8.2 Kpc, but such feeble differences 
may have gone undetected at the level of accuracy of our 
relative photometry (e.g. may be hidden in the "instrumen- 
tal width" of the observed sequences). 

6. The weak shoulder at [M/H] 0.2 in the MD of F2 is 

due to the handful of red stars around the ridge line of 
NGC 6553 discussed above, hence, at the present stage, 
cannot be trusted as a real feature of the MD (but it de- 
serves further investigation). 



3.4. Stellar Populations and radial distributions 

The stellar content of the outer region of M 33 is poorly ex- 
plored. It is interesting to investigate if the different kind of 
stars identified in our CMDs share the same spatial distribution. 
To check this point we defined three selection boxes that are de- 
picted in the Fl CMD shown in Fig. 7. Note that the boxes are 
defined in a range of magnitude in which incompleteness ef- 
fects should be weak or negligible (M/ < -3.0, e.g. more than 
2 mag above the limiting magnitude) and cover similar magni- 
tude ranges. 

As a guideline, we have superimposed three isochrones of 
solar metallicity and age = 20,40, 100 Myr (continuous lines) 
and three isochrones with [M/H] = -0.7 and age = 2 , 6, 12 
Gyr (empty symbols) from the set bv lGirardi etall d200(j) . The 
bluest box samples the upper MS, e.g. young stars with age 
< 100 Myr. The large box above the TRGB samples the bright 
AGB stars, the smaller box samples the brightest RGB stars. 
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Fig. 7. Description of the adopted selection boxes. The box 
around V — I = 0.0 selects a sample of young MS stars, the 
large box above the TRGB selects AGB stars, the small box 
below the TRGB selects RGB stars. The continuous lines are 
isochrones of solar metallicity, the o pen symbols are isoch rones 
at [M/H] = -0.7, from the set bv ldirardi etail ll2000h . The 
ages are reported in the plot. The small dots are the observed 
stars of Fl. 

Table 2. Distance Modulus, median Metallicity and standard 
deviations in different metallicity scales, and with different as- 
sumptions on the reddening. Note that the median and the mean 
metallicity are nearly coincident in all the considered cases. 





Fl F2 


Fl F2 


E(B-V) 


0.04 


0.08 


[Fe/H]zw 
o~ zw 

[Fe/H]co 

O~C0 

[M/H] 

C [M/H] 


-1.03 -1.03 
0.40 0.40 

-0.89 -0.89 
0.27 0.28 

-0.75 -0.74 
0.23 0.23 


-1.17 -1.18 
0.40 0.44 

-0.97 -0.98 
0.28 0.32 

-0.81 -0.81 
0.23 0.27 


(m - M)o 


24.64+0.15 


24.59±0.15 



Both kinds of tracers are associated with intermediate to old 
age stars, but are not necessarily linked thogether. In particular, 
brig ht RGB stars trace p opulations older than 1-2 gyr (see 
ISalaris. Cassisi & WeissL Eo02l and references therein). In the 
following, R must be intended as the projected angular distance 
from the center of the galaxy. 



In the upper panels of Fig . 8 the adopted selection boxes are 
superposed on the CMDs of Fl (left panel) and F2 (right panel). 
The lower panels of Fig. 8 show the cumulative radial distribu- 
tions of the stars falling in the boxes in the two different fields. 
In the radial range covered by Fl the distributions of RGB and 
AGB stars are indistinguishable. On the other hand, MS stars 
appear much more centrally concentrated since their distribu- 
tion see ms to end at R ~ 17' — 18', e.g. around 2 disc scale- 
lengths JVan den Berg hi ll 99 lh . According to a Kolmogorov- 
Smirnov test, the probability that the MS and RGB samples are 
drawn from the same parent population is P < 0.04%. This 
suggests that the RGB and AGB stars are not associated with 
the young disc component traced by MS stars. At least a (sig- 
nificant) fraction of them should belong to a more extended 
galactic component. 

The radial distribution of MS stars is not reported in the 
lower right panel of Fig. 8 since, in agreement with the above 
conclusion, F2 is virtually devoid of stars populating the MS 
box. It is surprising to note that in this field, AGB stars appear 
to follow a significantly different distribution with respect to 
RGB stars. The distribution of AGB stars is less centrally con- 
centrated and it is quite similar to uniform distribution on the 
sky. The latter fact would be naturally explained if the AGB 
sample of F2 would be dom inated by foreground contamina- 
tion (see also iDavidg e 120031). According to the predictions of 
the Galactic model by Robin et all d2003l) this seems to be the 
actual case. The model predicts that the number of Galactic 
stars falling in the AGB box is 48, less than compared to the 
61 actually observed. Hence ^ 80 % of the putative AGB stars 
in F2 are likely foreground stars. The impact is much smaller 
on the Fl stars where 584 stars are observed in the AGB box, 
hence the fraction of foreground contaminants is < 10 %. On 
the other hand the expected number of foreground stars falling 
in the RGB box is ~ 4, e.g. negligible in both fields. 

Therefore, the difference of radial distribution shown in the 
lower right panel of Fig. 8 is completely spurious. On the other 
hand, if we consider star counts and take into account the cor- 
rections for foreground contamination, it turns out that while 
the number of RGB stars drops by a factor 27 going from 
Fl to F2, the number of AGB stars decreases by a larger factor, 
42. This suggests that RGB stars may follow a more extended 
distribution with respect to AGB stars at large radii. 

All the above considerations seem to indicate that popula- 
tions of different characteristic ages follow different distribu- 
tions on the sampled scales, the older stars having more ex- 
tended distributions. This is suggestive of the presence of a 
weak "classical" old halo component in M 33. The only pre- 
vious indication in this sense (from field stars) is provided 
by the discove ry of a few candidate RR Lyrae variables by 
IPritchetl dl988h . while all other hints of the existence of an old 
spheroid al stellar component come from the study of g lobular 
clusters (Schomme ret all Il99 1 : Saraiedini et alll200 0). Note, 
however, that an extended and old disc component is also com- 
patible with our observations (see T04). 

In the upper panel of Fig. 9 we report the radial profiles 
of RGB and MS stars over the range(s) sampled by the present 
study. The MS profile is reasonably reproduced by an exponen- 
tial law with scale-length h = 9.2' up to R 17', in agreement 
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with the results bv lKentl ( ll987h . However, at R = 19' the sharp 
drop of the density already observed and discussed in Fig. 8 is 
clearly evident. On the other hand the RGB profile is well fit- 
ted by an exponential law with h = 4.9' over the whole radial 
range sampled by Fl. The density of RGB stars falls signifi- 
cantly below the adopted exponential profile in F2, suggesting 
a break in the observed profile in the range 20' % R ^ 25'. 
In our view, the most interesting result of this comparison is 
that, even ignoring the observed density cut-offs, the RGB and 
MS distributions do have significantly different profiles, again 
suggesting a different origin. 

The lower panel of Fig. 9 shows that in the radial range 
covered by Fl (10' < R < 20') the observed density profile 
of RGB stars is equally well fitted by the exponential law de- 
scribed above and by a R* law having effective radius R e = 2.7' 
(in fact, the R~ law provides a marginally better fit with respect 
to the exponential). It is interesting to note that the same R~ law 
was found to provide a good fit also to th e central bulge of M 33 
( Boul esteix et all Il980t iBothunl Il992l) . suggesting a possible 
connection be tween the bulge and the putative halo component 
(but see Step hens & FrogelH2002l for a detailed decomposition 
of the inner profile). While suggestive, the above result is lim- 
ited to the considered radial range (10' < R < 20') where the 
contribution of disc stars to the RGB population may be low. 
To correctly disentangle the contribution of the disc from that 
of the more extended component identified here, a complete 
sampling of the density profile from the center to the outskirts 
of the galaxy is needed, e.g. covering also the regions in which 
the surface brightness should be dominated by the exponential 
disc. This kind of analysis is clearly beyond the reach of the 
present study. 

4. Conclusions 

We have obtained a robust detection of the I magnitude of the 
TRGB in a field located at ~ 15' from the center of M 33, 
near the galaxy minor axis. Adopting the median metallicity 
we derived from the same data and the calibration of Mj RCB 
as a function o f the g lobal metallicity ([M/H]) provided by 
iBellazzini et alJ (|^)04), we have obtained a new estimate of 
the distance modulus of M 33, (m - M)o = 24.64 + 0.15. All 
the sources of uncertainty have been taken into account in the 
reported error bar. 

In Fig. 10, our distance modulus is compared with pre- 
vious estimates available in the literature. The large majority 
of the reported values are compatible, within the formal 1-cr 
errors, with our estimate. The only exceptions are provided 
by the oldest a nalysis of Cepheid variables, based on pho- 
togra phic plates ( Mad ore et all Il985t IChristian & Schommerl 
198 7b. If we exclude these two estimates as well as those 
bv lArgon etail Jl998l) and lGreenhill et all Jl993l) . that are af- 
fected by very large uncertainties, we obtain, from 19 inde- 
pendent estimates including ours, an average distance modulus 
(m-M)o = 24.69±0. 15 (average + standard deviation) in excel- 
lent agreement with our result. It is interesting to note that the 
standard error on the above average is just e — o~/ yl9 = 0.03 
mag. 




R [arcmin] R [arcmin] 

Fig. 8. Upper panels: the selection boxes described in Fig. 7 
are superposed on the CMDs of Fl (left) and F2 (right). The 
corresponding cumulative radial distributions are displayed in 
the lower panels. 



The photometric metallicity distributions described in 
Sect. 3.3 indicate that (a) the observed RGB population in the 
outskirts of M 33 has a typical metallicity ([Fe/H]zw — —1.0) 
that is intermediate between that of the halo of the Milky Way 
(\Fe/H] 7W ^ -1.5) an d that of M 31 ([Fe/H] zw * -0.6) (see 
IBellazzini et all 120031 for discussion and references); (b) the 
MDs are quite similar everywhere, within th e samp led regions, 
similar to the case of M 3 1 ( Bell azzini et alll2003l) . 

In the radial range 10' < R < 20', the AGB and RGB stars 
have a similar radial distribution, much more extended than that 
of young MS stars which appear to decrease abruptly around 
R - 18'. The distribution of RGB stars is equally well fitted 
by an exponential law (but not the same that fits the distribu- 
tion o f MS stars) or by the sa me R 1 ^ 4 law that fits the central 
bulge (Boulesteix et al., 1980). The density of RGB stars is ob- 
served to fall far below that predicted by both the adopted best- 
fit models at R ~ 27', but the actual break may occur anywhere 
between R = 20' and R = 27' (e.g. between 4 and ~ 5 expo- 
nential scalelen ghts, in good ag r eemen t with the complete pro- 
file obtained by Fer guson et all (12 004) from the huge database 
presented in iMcConnachie et all 1120041) : A. Ferguson, private 
communication). The above results may suggest that a weak 
old-halo component is indeed present also in M 33, but proba- 
bly the point could be established only with an analysis of the 
stellar kinematics in the considered radial ranges. The possibil- 
ity to interpret the observations of the stellar content at large 
distance from the center of the galaxy in terms of an unexpect- 
edly metal rich halo or of a very extended old disc is (again) 
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Fig. 9. Upper panel: stellar density profiles for MS (open cir- 
cles) and RGB stars (open squares) compared with two dif- 
ferent exponential laws. All the estimates with R < 20' are 
computed over sections of concentric annuli 2' wide. The RGB 
points at R = 27.6' is the density estimate obtained from F2 as 
a whole. Lower panel: expansion of the RGB density profile. 
The continuous line is the same exponential law displayed in 
the upper panel, the dotted line is a R+ law with R e - 9.2'. The 
associated uncertainty is the Poisson noise of the star counts 
propagated to the adopted unities. 

reminescent of the case of M 31 (see iFerguson & Johnson! 
l2QQlUBellazzini et allEooi for discussions and references). 
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